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Abstract. The escape probability is a deterministic tool that quantifies some aspects 
of stochastic dynamics. The goal of the present work is to analytically examine escape 
probabilities for dynamical systems driven by symmetric a-stable Levy motions. Since 
escape probabilities arc solutions of a type of intcgro-diffcrcntial equations (i.e., differen- 
tial equations with nonlocal interactions) , asymptotic methods arc offered to solve these 
equations to obtain escape probabilities when noise is sufficiently small. Two examples 
are presented to illustrate the asymptotic methods, and asymptotic escape probability 
is compared with numerical simulations. 



1. Introduction 

Stochastic dynamical systems arise as mathematical models for complex phenomena in 
biological, geophysical, physical and chemical sciences, under random fluctuations. Unlike 
the situation for deterministic dynamical systems, an orbit of a stochastic system could 
vary wildly from one sample to another. It is thus desirable to have efficient tools to 
quantify stochastic dynamical behaviors. The escape probability is such a tool. 

Non-Gaussian random fluctuations are widely observed in various areas such as physics, 
biology, seismology, electrical engineering and finance [T6|[9|[ll]. Levy motions are a class 
of non-Gaussian processes whose sample paths are discontinuous in time. For a dynamical 
system driven by Levy motions, almost all orbits are discontinuous in time. In fact, these 
orbits are cadlag (right-continuous with left limit at each time instant), i.e., each of these 
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orbits has countable jumps in time. Due to these jumps, an orbit could escape an open 
domain without passing through its boundary. In this case, the escape probability is the 
likelihood that an orbit, starting inside an open domain D, exits this domain first by 
landing in a target domain U in D c (the complement of domain D). 

For brevity, in this paper we only consider scalar stochastic dynamical systems. Let 
{X t ,t ^ 0} be a real-valued Markov process defined on a complete filtered probability 
space J 7 , {J-'tjtzoi P)- Let D be an open domain in R. Define the exit time 

t D c := inf{t > : X t E D c }, 

where D c is the complement of D in M. Namely, t^c is the first time when X t hits D c . 

When X t has cadlag paths which have countable jumps in time, the first hitting of D c 
may occur either on the boundary 3D or somewhere in D c . For this reason, we take a 
subset U of the closed set D c , and define the likelihood that X t exits firstly from D by 
landing in the target set U as the escape probability from D to U, denoted by p(x). That 
is, 

p{x) = ¥{X TDC G U}. 

If Xt is a solution process of a dynamical system driven by a symmetric a-stable Levy 
motion, by [10J [7J, the escape probability p(x) solves the following Balayage-Dirichlet 
"exterior" value problem: 

Ap = 0, xeD, 
p\dc = f, 

where A is the infinitesimal generator of X t and tp is defined as follows 

, . _ / 1, xeU, 
v[x) ~ \ 0, x e D c \ U. 

However, Eq.((T]) is usually an integro-differential equation and it is hard to obtain exact 
representations for its solutions. Here we use asymptotic methods to examine its solutions. 
More precisely, (i) for a dynamical system driven by a Brownian motion combined with a 
symmetric a-stable Levy motion, an asymptotic solution of Eq.((T]), or escape probability 
p(x) from D to U, is obtained by a regular perturbation method; (ii) for a dynamical 
system driven by a symmetric a-stable Levy motion alone, the escape probability p(x) is 
obtained by a singular perturbation method. 

This paper is arranged as follows. In Section [2j we introduce symmetric a-stable Levy 
motions and their infinitesimal generators. In Section [31 a regular perturbation method is 
applied to examine escape probability for dynamical systems driven jointly by Brownian 
motion and symmetric a-stable Levy motions. Escape probabilities for dynamical systems 
driven by symmetric a-stable Levy motions alone are studied in Section H] by a singular 
perturbation method. Two examples are presented in Section In the Appendix (Section 
E]), solvability of a type of integro-differential equations is discussed. 
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2. Preliminaries 

In this section, we recall basic concepts and results that will be needed throughout the 
paper. 

Definition 2.1. A process Lt — (^t)t^o with Lq — a.s. is a Levy process or Levy motion 
if 

(i) L t has independent increments; that independent of L v — L u if (u, v) D 
(M) = 0; 

(ii) L t has stationary increments; that is, L t — L s has the same distribution as L v — L u 
if t — s = v — u > 0; 

(Hi) L t is stochastically continuous; and 
(iv) L t is right continuous with left limit. 

The characteristic function for L t is 

E (exp{izL t }) = exp{t#0)}, zGl. 

We only consider scalar Levy motions here. The function \1/ : R — > C is called the 
characteristic exponent of the Levy process L t . By Levy-Khintchine formula, there exist 
a nonnegative number Q, a measure v on K satisfying 

z/({0}) = 0and / (|w| 2 A l)u(du) < oo, 

JR\{0} 

and also a real number 7 such that 

*(z) = t 1 z-\Qz 2 +f (e^ tt >-l-i(z,«>l M<1 )i/(d«). (2) 

1 JR\{0} 

The measure v is called a Levy jump measure, Q is the diffusion, and 7 is the drift. 

We now introduce a special class of Levy motions, i.e., the symmetric a-stable Levy 
motions Lf . 

Definition 2.2. For a G (0,2). A scalar symmetric a-stable Levy motion Lf is a Levy 
process with characteristic exponent 

q( z ) = -\ z \ a , zeR. 

Thus, for a scalar symmetric a-stable Levy motion L", the diffusion Q = 0, the drift 
7 = 0, and the Levy jump measure v is given by 

u(du) = - — -r7-dw, 

V I \ u \l+a 

where 

ar((l + q)/2) 
1,Q 2 1 — 7r 1 /2r(l -a/2)' 
Let Co(M) be the space of continuous functions / on M satisfying lim f(x) — with 

\x\—>oo 

norm ||/||c (R) = sup|/(x)|. Let Cq(1R) be the set of / G Co(M) such that / is twice 
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differential) le and the first and second order derivatives of / belong to Co(M). Let C a be 
the infinitesimal generator of . By [131 Theorem 31.5], 



(C a f)(x) = / (f(x + u)- f(x) - f'{x)ul luHl )is(du), 

JR\{0} 

where / G C 2 (R). For any e > 0, eL" is also a scalar symmetric a-stable Levy motion, 
and its Levy measure is e (B) = v{-B) for B G <S^(R) (Borel a-algebra on R). Thus, its 
infinitesimal generator is 

(C e J)(x):= [ (f(x + u)-f(x)-f'(x)ul lu]<i y(du), /GC 2 (R). 

</R\{0} 

Applying the representation of is, one can obtain that 

(C e J)(x) =e a [ (f(x + u)- f(x) - f'(x)ul luHl )is(du). 



3. Escape probability of a SDE with Brownian motion and symmetric 

a-STABLE Levy motions 

Let {W(t)} t ^>o be a scalar standard ^-adapted Brownian motion, and Lf a scalar 
symmetric a-stable Levy motion with a G (0, 2) and independent of W t . Consider the 
following scalar stochastic differential equation, with drift b, diffusion a and intensity 
e(> 0) of Levy noise, 

dX t = b(X t ) dt + a(X t ) dW t + edL«, . 
X = x. W 

Assume that the drift b and the diffusion 0) satisfy the following conditions: 
(Hfe) there exists a constant Cf> > such that for x, y G R 

\b(x) - b(y)\ sC C b \x - y\ ■ log(\x - y\~ l + e); 

(H CT ) there exists a constant C a > such that for x, y G R 

- a(y)| 2 ^ C CT |x - y| 2 • logflx - y\~ l + e). 

Under (H;,) and (Ho-), it is well known that there exists a unique strong solution to 
Eq.(j3])(see [12] )■ This solution will be denoted by X t (x). By Theorem 3.3 in [7J, the 
escape probability p{x) for X t (x), from D = (A, B) to C/ = [5, oo), satisfies the following 
integro-differential equation 

b(x)p'(x) + ^cr 2 (x)p"(x) + e a I (p(x + u) — p(x) - I\ u \^ip' (x)u) u(du) = 0, 

2 JR\{0} 

xe(A,B), (4) 

and the 'exterior' conditions 

P(x)\(-cc,A] = °> P(^)|[B,oo) = L (5) 
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We consider the solution for p(x), when e > is sufficiently small. Assume that p(x) 
has the following regular expansion 

p(x) = p (x) + e a pi(x) + e 2a p 2 (x) H . (6) 

Substituting (JS)) into @ and equating like powers of s, we obtain a system of equations 
for the recursive determination of Pj(x). The leading-order equation for po(x) is 

b(x)p' (x) + -a 2 {x)p' Q {x) = 0, x G (A, B) (7) 
with boundary conditions 

p (A) = 0, po(B) = 1. (8) 
Using the boundary conditions, we solve Problem (]7j) and (jHJ) to get 



p [x) 

where 4>(u) := 2b(u) / a 2 (u) . 
Next, the equation for p\{x) is 

«.)ri<.) + + / (»<* + «) - »<*) - WiM-) ^du) = 0, 

1 </R\{0} V 7 

xe(A,B), (9) 

with boundary conditions 

Pi(A) = 0, Pi(B) = l. (10) 



9( x)==/ (,( I+ »)- a W-/,**W«)*). 



Set 



'M\{0} 

Then Eq.Q is transformed into the following equation 

b(x)p' 1 (x) + -a 2 (x)pl(x)+g(x) = 0, xe(A,B). (11) 



By solving Problem fllip and (11 01) we get 



A V 

-p [x) I e~ ^<t>( u ) du ( J — eJA^'^au ) as + p {x) 



Thus we have an asymptotic expression for escape probability, i.e., solution of Eq.fllJ), 
for e sufficiently small, 

p(x) fa p (x) + e a pi(x). (12) 
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4. Escape probability of a SDE with symmetric ^-stable Levy motions 

Consider the following stochastic differential equation with a symmetric a-stable Levy 
motion, with 1 < a < 2, on K 

' dX t = b(X t ) dt + edLf , (13) 

where the drift b satisfies (H&). 

By [T5l Theorem 3.1], Eq. ffTB"]) has a unique solution X t (x). From Theorem 3.3 in [7], 
the escape probability p(x), for X t (x) from D = (A,B) to U = [B,oo), satisfies the 
following integro-differential equation 



b(x)p'(x) + e a / (p(x + it) — p(x) — p'(x)u)u(du) = 0, 

JR\{0} 

xe(A,B), (14) 

with the 'exterior' conditions 

P(a0l(-«v4] = 0, (15) 

p(x)\[B l0 o) = 1. (16) 

We now try to construct an asymptotic solution of (fl4"]) . (|15l) and (fl6l) for sufficiently 
small s > 0. We consider the following four different cases, depending on the dynamical 
behavior of the corresponding deterministic dynamical system x = b(x). 

Case 1: b(x) > for x G (A, B). In this case the deterministic dynamical system 
x = b(x) has no equilibrium states and all orbits move to the right. 

Thus it is reasonable to require that p(x) — > 1 as e — > 0. We assume that p(x) has the 
following expansion 

p(x) = p (x) + e a pi{x) + e 2a p 2 (x) + ■■■ . (17) 

Substituting (fTTj) into ( IT4|) and equating like powers of e, we obtain a system of equations 
for the recursive determination of Pj(x). The leading-order equation for po(x) is 

b(x)p' (x) = 0, SG(AB), 

and thus Po( a; ) = 1 for a; G (A, £>), because p(rr) — > 1 as e — > 0. Since po^) does not satisfy 
the boundary condition (IT5|) . it is necessary to construct a boundary layer correction to 
Po{x) near a; = A. 

We introduce a stretched variable 

x- A 

f = — 

with /3 > determined later. Defining -F(^) = Po(A + ^'e' 3 ) and inserting it into Eq.( IT4l) . 
we obtain 

b(A + &P)F'(£)e-P + e a - a P [ (F(£ + u) - F(£) - F'(t)u)v(du) = 0. (18) 



Set — j3 = a — a(3. That is, we take j3 = -^-r. Multiplying Eq.f ll8j) with e 13 and letting 
e — > 0, we imply that 

b(A)F'(0 + / + u) - F(0 - F'({)u)v(du) = 0, (19) 

Jr\{o} 

with the boundary condition 

F(€)=0, £<0, (20) 

and the matching condition 

lim F(0 = 1. (21) 

As seen in Appendix at the end of this paper, we know that the system (fT9~]) - (j2T|) is 
solvable, although the solution cannot be expressed in terms of elementary functions. So, 

. . ^ fx - A s 
Po{x) = F — — 

Thus an asymptotic solution of p(x) is, for sufficiently small e, 

' x - A" 



p{x) 

£a-l 

Case 2: b(x) < for x G (A, B). Again, in this case the deterministic dynamical 
system x = b(x) has no equilibrium states and all orbits move to the left. 
Thus as e 0, p(x) —> 0. We assume that p(x) has the following expansion 

p(x) = p (x) + e a Pl (x) + e 2a p 2 (x) + ■■■ . (22) 

Similar to Case 1, we obtain the leading-order equation for po(x) 

b(x)p' (x) = 0, i£ (A,B). 

So, po(x) = for x G (A, B), because p(x) — > as e — > 0. Since po(x) does not satisfy 
the boundary condition ffT6]) . it is necessary to construct a boundary layer correction to 
Po(x) near x = B. 

We introduce a stretched variable 

B — x 

where (3 is the same as one in Case 1. Defining G(s) = Po(B — qe 13 ) and inserting it into 
Eq.( 1T4"j) . we obtain 

- b(B - qe^GX^e-P + e a ~ afi [ [G(q - u) - G{q) - G'^)(-u)]u(du) = 0. (23) 

</R\{0} 

Multiplying Eq.(J23J with and letting e — V 0, we obtain 

- b{B)G'{q) + f (G(q + u) - Gfc) - G'(s)u)u(du) = 0, (24) 
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with the boundary condition 

G{q) = l, ^0, (25) 

and the matching condition 

lim G(q) = 0. (26) 

As seen in Appendix, the system fl24l) - fl26l) is solvable, although the solution cannot be 
expressed in terms of elementary functions. So, 

p (x) = G ( — a X 

Thus we obtain an asymptotic solution of Eq. (1141) 

B-x 



p(x) ~ G 

Case 3: There exists a x G (A,B) such that b(x) = and b'(x) > 0. In this case the 
deterministic dynamical system x = b(x) has one unstable equilibrium state x. Then as 
e 0, p(x) — > 1 for x < x ^ B and p(x) — > for A ^ x < x. We assume that p(x) has 
the following expansion 

p{x) = p (x) + e>i(x) + e 2a p 2 (x) + ■■■ . (27) 

As in Case 1, we obtain that the leading-order equation for po{x) is 

b(x)p' (x) = 0, x E (A,B). 

So, 

1, x < x ^ B, 
0, A ^ x < x. 

Although po(x) partially satisfies the 'exterior' conditions (TTBT) and (1761) , the value of po( a; ) 
around x is unknown. Therefore, it is necessary to construct an internal boundary layer 
correction to po(x) near x = x. 
We introduce a stretched variable 

x — x 



p (x) 



V = • 

e 

Define H(rj) = po(x + rje) and insert it into Eq. (fl4j) . Then Eq. fll4l) becomes 

b(x + ^H'i^e- 1 + / (H(r) + u) -H{rj) - H'(r))u)v(du) = 0. 
Jm\{o} 

Letting £ — > and using the L'Hospital's rule, we get 

b'{x)r)H'(r})+ I (H{r) + u) - H(r}) - H'(r})u)u(du) = 0, (28) 

JR\{0} 

with the matching conditions 

lim H(r]) = 0, (29) 

•q— y— oo 



lim H(rj) = 1. (30) 

77— >oo 



As seen in Appendix, Eq.( l2"gj) is solvable. So, 

x — x 

e 

Thus we obtain an asymptotic solution of Eq. ([T 



Po(x) = H 



p(x) ~ if 



x — X 



Case 4: There exists a x G (A, 5) such that b(x) = and b'(x) < 0. In this case the 
deterministic dynamical system x = b(x) has one stable equilibrium state x. We further 
require that b(A)b(B) 7^ 0, i.e., A and B are not equilibrium states for x = b(x). 

We assume that p(x) has the following expansion 

p(x) = po(x) + e a jn(x) + e 2a p 2 (x) + ■■■ . (31) 

As in Case 1, we obtain the leading-order equation for po(^) 

b(x)p' (x) = 0, x E (A,B). 

So, po{x) = C for x G (A, B). Because of not knowing at which endpoint there will be a 
boundary layer correction, we construct asymptotic approximations near both endpoints. 
If there is a boundary layer correction near x = A and x = B, respectively, as in Case 1 
and Case 2, we get near x = A 

x - A 



B-x 



Po(x) = C-F 

and near x = B 

p (x) = C + (1 - C) ■ G 

Thus we have an asymptotic solution 

p(x) ^C-F + (1 - C) ■ G 

Since F, G cannot be expressed in terms of elementary functions, it is hard to determine 
C. But for a concrete example in the next section, Example I5.2[ we introduce a method 
to determine the value of C. 



5. Examples 

In this section we consider two examples. Example 15.11 and Example 15.21 correspond to 
our methods in Section E] and Section HJ respectively. 
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Example 5.1. Consider the following scalar SDE with a Brownian motion and a sym- 
metric a-stable Levy motion: 

' dX t = dW t + edL?, 
X = x. 

The unique solution is denoted as X t (x). We take (A, B) = (—1, 1) and [B, oo) = [1, oo). 
The escape probability p{x), for X t (x) from (—1,1) to [l,oo), satisfies the following 
integro- differential equation 



lp"(x)+e a 



(p(x + u) -p(x) - I\ u \ <x p' \x)u) i/(du) = 0, 

xe (-1,1), 



(32) 



and the exterior conditions 

p(s)|(_ 00> _i] = 0, p(x)|[i i0o) = 1. 
We seek an asymptotic solution of p(x) as follows 

p(x) m po(x) + e a p 1 (x), 

where 



and 



Pi(x) 



Cl a 



p Q (x) 



a)(3-a) [(-*■ 



2±i, -1<X<1, 

0, x ^ -1, 

1, x^l, 



X 



,3-a 



(-a)(l-a)(2-a)(3-a) 
x±l Cl 

2 (-a)(2-a)(3-a)^ ^ 2 ' 1 ^ X ^ A 



2 3-a + (3 _ a )2 2 ^ Q (x + !)-(! + x) 3 ~ a l 



0. 
1. 



X -1, 
X > 1. 



iVeatf we use a numerical method to study Eg. / fffjgj) . i??/ t/je exterior conditions, Eq. 
is changed into 



-p (x) - e a C ha p (x) - e a C^ 

Z a — I 



(l + x)~ a + (l-x)~ a 



-e a C h 



1 p{x + u) — p(x) — p (x)u 



a a 

x+l 



p(x) 



x-l 



it 



1+a 



du + e a d a 



p(x + u) — p{x) 



\ll 



1+a 



du 



+e a C ha {1 X) ~° = 0, 



a 



/or — 1 < x < and 



\p"{x) - e^d,, 1 (] , X) ' a+1 p'(x) - e a C 1>a 
2 1 — a 



(l + x)- Q 



a 



p(x) 



-i 



p(x + u) — p(x) 



u 



du + e Q C 1 



p(x + u) — p(x) — p \x)u 



\u 



1+a 



du 
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+e°Ci J 1 ^ = 0, 
a 

for ^ x < 1. For any fT G N, dmii equally the interval [—2,2] into subintervals 
and let x k := fc/i and P k := p(xfc) /or — 2K ^ k ^ 2K, where h — j?. Thus, by central 
difference for derivatives and punched-hole trapezoidal rule, one can obtain 

1 iVi - 2P k + P k+l l-tl + x^ 1 ^-^ 

2 /i 2 £ 1,Q a - 1 2/i 



a a J a 

P i+k - P k - {P k+1 - P k -i)xi/2h 



+e a C ha hY 

^— ' i=-K~k \Xi\ 1+a 

t£ -wr"^ = o, 



/or fc = — K + !,••■ , — 1 and 



1 Pfc-i - 2P fc + P fc+1 1 - (1 - a: fc )~ a+1 P fc+1 - Pfc-i 

2 /i 2 E 1,Q 1 - a 2/i 



a a a 



+6 C 1>a h^ i= _ R ^ 

+^c 1Q / i v' _/ ' Pm ,r/ fc =o, 

^— ' i=-K-k \Xi\ l+a 

for k = 0, 1, • • ■ , K — 1. Fere £/ie modified summation symbol "^2" means that the 
quantities corresponding to the two end summation indices are multiplied by 1/2 waen 
£/iey are added. By the error analysis in [3], the following scheme have 2nd-order accuracy 
for any < a < 2 

P fe _ x - 2P fc + P k+l 1 - (1 + x k )~ a+1 P k+1 - P fc _ a 

Cfc /r 2 " 5 Gl ' a cT^l 2h 

^ r f (l+x k )- a , (l-x k )- a \ p (1 - 

— £ ^l,a 1 Pfc + £ ^1,< 



a a J a 

, «^ Pi+k ~ Pk ~ {Pk+i ~ Pk-i)xi/2h 

-+E Oi a Ai > ; — 

^— ' i=-K-k \Xi\ 1+a 

,__/K-fc P. , _ P, 

%# = 0, (33) 



for k = —K + !,••• , — 1 and 



P fc _! - 2P fc + P fc+1 1 - (1 - X fc )- Q+1 P, +1 - P fc _! 

h 2 1 — a 2h 
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(a) (b) 




Figure 1 . Comparison between the asymptotic solution and the numerical 
solution of Eq. (l32|) for small e. (a) a = 0.5, e = 0.01. The asymptotic 
solution is shown with dashed line while the numerical solution is displayed 
with solid line, (b) a = 1.5, e = 0.01. (c) a = 0.5, e = 0.001. (d) a = 
1.5, e = 0.001. 



■ Xk )- a (i-x k r a \ (i-x k ) 

-E ^1,q 1 -ffe + £ W,' 



a a J a 



+e C ha h^ i= _ k 

+ £ «C 1 a fcV P ; +fc ,~ Pfc = 0, (34) 

' i=-K-k \Xi\ 1+a 

for k = 0, 1, ■ • • , K — 1, where Ch = | — £ Q Ci iQ C(a — l)h 2 ~ a and ( is the Riemann 
zeta function. By solving these linear equations (CHJ) and (34\ ), we obtain the numerical 
solution of Eq.ftM}. 

Figure U\ shows that if e is fixed and a turns large, the difference between the asymp- 
totic solution and the numerical solution of Eq. ( TffDj) will become small; if a is fixed and e 
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becomes large, the difference will turn large, because the asymptotic solution is for suffi- 
ciently small e. 

Example 5.2. Consider the following scalar SDE with a symmetric a-stable Levy mo- 
tion, with 1 < a < 2, 

dX t = -X t dt + edLf, 
X = x. 

The unique solution is denoted as X t (x). We take (A,B) = (—1, 1) and [B, oo) = [1, oo). 
The escape probability p(x), for X t (x) from (—1,1) to [l,oo), satisfies the following 
integro- differential equation 



— xp'(x) + e a / (p(x + u) — p(x) — p'(x)u) z/(du) = 0, 
Jr\{o} 

xe(-l,l), (35) 

with the exterior conditions 

p(s)|(_ 00> _i] = 0, p(x)|[i i0o) = 1. (36) 

Since b(x) = -a:, 6(0) = 0,1/(0) < 0,6(-l)6(l) = -1 ^ 0, by the result of Case 4 in 
Section^ an asymptotic solution of p(x) is given by 

p(x) ^C-F (^-) + (1 - C) ■ G (^-) . (37) 



Specially, take the Levy measure 

v ( du ) = ■ liui^du, 

where k > is a constant (\16\). Thus, the function F can be given explicitly by 

F(x) 

where 7 > satisfies the following integral equation 

7 - f (e-^ - 1 - (-j)u) -J^du = 0. 



1 - e~^ x , x > 0, 
0, x < 0, 



By the relation between F and G, we can obtain 



G(x) = I I 



t x , x > 0, 
x < 0. 



So, the asymptotic solution of Eq. i fffS)) is given by 

p{x) « C (l - exp (-7 + C) exp (-7 / 
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To determine C , we multiply Eq. by the solution p(x) of the steady Fokker-Planck 
equation 

- (-xp(x)Y + e a I (p(x + u)- p(x) - p'(x)u) u{du) = (38) 

JR\{0} 

and integrate over (— 1, 1), to obtain 

( — xp{x))p' {x)dx + e a I p{x)dx / (p(x + u) — p{x) — p'{x)u) v{du) = 0. (39) 



-l J -i 

To by integration by parts and using ( ESP ? we 9 e ^ 

—p{l)—e a I p(x)dx / (p(x + u) — p{x) — p'(x)u) v(du) 

J-l JR\{0} 

+e a / p(x)dx / (p(x + u) — p(x) — p'(x)u) u(du) = 0. 

J -l </R\{0} 

Applying Cauchy principal value and IfSh}) . we have 



1 



p(x)dx 



(40) 



-B?/ [D Proposition 3.2], £/ie Fourier transform of p(k) is given by 

p(k) = exp{--\k\ a }. 

a 

Replacing p(x) and p(x) by ^- f R e zxk p(k)dk and (3l\ ), respectively, and letting e — )■ 0, we 
can obtain C from 

6. Appendix 

In this appendix, we consider the system (|T9l)-(l20l)-( l2T|) . arising in the analysis in this 
paper, including a useful identity, solvability condition, and the matching condition for 
the asymptotic solution. 

6.1. A useful integral identity. First, we prove a useful integral identity: 

- 1 " (-A)«) ~^du = A 7 2 7 } . (41) 
Proof. By integration by parts, 

(e-*« - 1 - (-X)u) -^du = -- / (e-^ - 1 - (-A)u) d«- 
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— / u- a (e~ Xu - l)du 
a Jo 



A 2 


a\ 


[a- 


1) 




X a 




a\ 


[a- 


1) 




\ a 




a\ 


'a - 


1) 



oo 

e -Xu u -a+l du 

oo 

e- u u~ a+l du 



r(2-a). 

This completes the proof. □ 

6.2. Solvability. Next, we consider solvability of the system fTr^ - (l2T)]) - fl2T]) . By jH El 
E], we know that Eg. f TT9|) is solvable, but the solution cannot be expressed in terms of 
elementary functions. 

6.3. Asymptotic behavior. Finally, we examine the asymptotic behavior of the solution 
to the system (fT9l)-(|20l)-(|2T|). 

For £ > 0, using ( |20|) . Eq.(fl9l) is transformed to 

b(A)F'(C) + / * - F'(On) u(du) 

J — oo 

/oo 
^ {F{i + u)- F(0 - F'(0«) u(du) = 0. (42) 

Denote the Laplace transformation of F(£) by -F(A), i.e. 



F(A) 



POO 

/ e- A «F(0d^, A>0. 
Jo 



Now applying the Laplace transformation to Eq. (1421) . we get that for the first term of its 
left hand side (LHS) 

/*oo /*oo 

/ e- x ^b(A)F'(0^ = b(A)X / e" A «F(0d£ 
Jo Jo 

= b{A)XF(X). (43) 

For the second term in the LHS of Eq. 

r>00 / /■— £ 



J™ e" A « (| (-F(C) - u(du)j d£ 

/■oo /•-§ ,>oo 

- / e^ A «F(0de / u{du) - / e- A «F'(0^ / uv{du) 

JO J —oo JO J —oo 



-oo JO J —oo 



poo /*°° r~ ^ ^ 
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POO f~i 

+ / e-*F(0^d*. (44) 
Jo s 

Now we deal with the third term in the LHS of Eq. (|42p . Firstly, let £ + u = u and then 

/oo 
^ (F(^ + u)-F(0-F'(0«)Kd«) 

= / (m - F(0 - F'(0(u - 0) , _ % +a du. (45) 
./o I" si 

Secondly, applying the Laplace transformation to Eq. (|45l) . we have 



oo 



where 



=: h-h-h, (46) 



POO poo S~1 



\u-tr 

poo poo s~i 

'■=/ ^ d? i F( «^ ds - 

/•oo /*oo r < 

For /i, interchanging the order of integrals, we obtain 



/*00 /'OO 

h = F(u)du / e 

■/o Jo 



oo 

e- Xu F(u)du I e~ A(g -") d£ 



n- s ^| 1 + Q 

C*l,a 

T+a ' 



/oo /*oo 
e" A "F(M)dM J e- Xu j^^du 

/0 /~i 
H 1+a 

r*oo noo fi 

+ / e~ M F(u)du / e~ Xu f^du. (47) 



\u 



We now treat /2 as follows: 



/ 2 = / e-*F®dZl -%-d« 



o |«-£| 1+a 

/*oo poo r i 

/OO 
-« M 1+Q 
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oo 



>r^F(u)du I ° ^du 



\u 



"OO pOO f i 

-I e~^F(u)du r^du. 



For h, we have 



du 



roc roc 
poo poo f < 



p OO p OO y-Hf 

poo /Hf 

+ 1 e-*F(0^de, 



where the Cauchy principal value is used in the last equality. 

Substituting (jl7|). (H5 jl and p?]) into (gSJ), and using (JH]), we obtain 



OO 



/^t /*oo /*oo A v 

r^F(u)dul e^^du + e-^F(u)du e^-^c 



r*0 s~i poo poo /~i 

^F(u)dul ^du- e-^F[u)du ±£d* 



/oo poo S~1 



OO /"OO 

e 



oo rO 



^F(u)du I (e" A « - 1 - (-A)u) ^du 







\U 



+ / e-^F(n)du / ( e ^« - 1 - (-A)u) ^du 

poo poo r < poo 

-2jf e-*F®dtj Xu^du-J^ e-*F(0^d£ 

poo pO r < 

\ e-^F(u)du / (e" A « - 1 - (-A)u) ^du 

JO J-U \ U \ 

/oo poo r < poo /~i 

e-*F(m] ^^p-J-l 



a{a — I) 
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Combining (H3]) ; (EH]) and ([50]), we get 

b(A)XF(\) - / e"*F(Od£ / f^du - / e -*F(£)d£ / Au^du 

JO J-oo \ u \ JO J£ \ u \ 

+ [ e- X *F(u)du £ (e^« - 1 - (-A)«) ^du + ^."^ ^(A) 



e-^F(u)du [° (e" A « - 1 - (-A)u) ^du 



+ / e-*F(£)d£ / - 1 - A«) ^d M 
Jo Jo Pi 



= 0. 
Therefore, 



^ (A) a(a - l)b(A)X + A«r(2 - a)C ha ' ' 6 

* e Au - 1 - Am , 
du 



1 + Aw, 
-du 



l+a 



It 



l+a 



Thus as A — > 0, -F(A) — » oo, which implies that F(£) approaches a constant for £ 3> 1 and 
satisfies the matching condition (|21[) . 
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